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基础 研究 
基于 投影 数据 全 广义 变 分 最 小 化 的 低 剂 量 CT 重 建 


牛 善 洲 , 吴 恒 , 喻 泽 峰 , 郑 子 君 , 喻 高 航 
闽南 师范 大 学 数学 与 计算 机 科学 学 院 , 江 西 赣州 341000 


摘要 :目的 提出 基于 投影 数据 全 广义 变 分 最 小 化 的 低 剂量 CT 重建 方法 。 方 法 首先 ,通过 非 线 性 Anscombe 变换 将 满足 
Poisson 分 布 的 投影 数据 转化 为 近似 Gaussian 分 布 ,然后 基于 全 广义 变 分 正则 化 模型 对 变换 后 的 Gaussian 型 数据 进行 噪声 去 
除 。 最 后 ,对 去 噪 的 结果 进行 Anscombe 逆 变换 后 实现 传统 的 滤波 反 投影 (EBP)CT 重 建 。 结 果 数值 体 膜 实验 结果 表明 本 文 提 
出 的 方法 可 以 大 大 地 改进 重建 图 像 的 质量 。FBP 方 法 重建 的 Clock 和 Shepp-Logan 体 膜 图 像 的 信 噪 比分 别 为 17.752 dB 和 
19.379 dB ,本 文 方法 重建 的 图 像 的 信 噪 比 提 高 到 24.0352 dB 和 23.4181 dB, FBP 方 法 重建 方法 重建 的 Clock 和 Shepp-Logan 
体 膜 图像 的 均 方 误差 分 别 为 0.86% 和 0.58% ,本 文 方法 重建 的 图 像 的 均 方 误差 降低 到 到 0.2% 和 0.23%。 结 论 本 文 方法 可 以 在 
影 数据 不 满足 分 段 常数 假设 的 前 提 下 去 除 噪 声 和 条 形 伪 影 ,从 而 提高 低 剂量 CT 图 像 重建 质量 。 
关键 词 : 低 剂量 CT 重建 ;全 变 分 ;全 广义 变 分 正则 化 ;Gaussian 分 布 ;滤波 反 投影 算法 
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Total generalized variation minimization based on projection data for low-dose CT 
reconstruction 


NIU Shanzhou, WU Heng, YU Zefeng, ZHENG Zijun, YU Gaohang 
School of Mathematics and Computer Science, Gannan Normal University, Ganzhou 341000, China 


Abstract: Objective To obtain high-quality low-dose CT images using total generalized variation regularization based on the 
projection data for low-dose CT reconstruction. Methods The projection data of the CT images were transformed from Poisson 
distribution to Gaussian distribution using the linear Anscombe transform. The transformed data were then restored by an 
efficient total generalized variation minimization algorithm. Reconstruction was finally achieved by inverse Anscombe 
transform and filtered back projection (FBP) method. Results The image quality of low-dose CT was greatly improved by the 
proposed algorithm in both Clock and Shepp-Logan phantoms. The signal-to-noise ratios (SNRs) of the Clock and Shepp- 
Logan images reconstructed by FBP algorithm were 17.752 dB and 19.379 dB, which were increased by the proposed algorithm 
to 24.0352 and 23.4181 dB, respectively. The NMSE of the Clock and Shepp-Logan images reconstructed by FBP algorithm was 
0.8676 and 0.5876, which was reduced by the proposed algorithm to 0.276 and 0.2376, respectively. Conclusion The proposed 
method can effectively suppress noise and strip artifacts in low-dose CT images when piecewise constant assumption is not 
possible. 

Keywords: low-dose CT reconstruction; total variation; total generalized variation; Gaussian distribution; filtered back-projec- 
tion algorithm 
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辨 率 上 的 日 越 表 现 , 已 经 广泛 应 用 于 临床 诊断 和 治疗 。 


低 剂量 的 CT 扫描 方案 虽然 能 在 一 定 程度 上 减少 X 射 线 
的 照射 剂量 ,但 采集 到 的 投影 数据 会 含有 较 高 的 噪声 ， 


CT 图 像 的 质量 与 X 射 线 辐射 剂量 密切 相关 ,剂量 越 高 
图 像 质 量 越 好 ,然而 ,过 量 的 X 射 线 照射 又 会 诱发 恶性 
肿瘤 .白血病 以 及 其 他 遗传 性 疾病 和。 当前 ,临床 上 采用 


直接 滤波 反 投 影 重建 图 像 的 质量 往往 无 法 满足 临床 诊 
断 的 需要 ”。 因 此 ,如 何在 降低 X 射 线 辐射 剂量 的 同时 
保证 图 像 质 量 已 经 成 为 CT 成 像 领域 的 迫切 要 求 。 

为 了 降低 X 射 线 的 辐射 剂量 ,许多 优化 的 低 剂 量 
CT 扫描 方案 ”以 及 抑制 噪声 和 伪 影 的 重建 算法 “” 相 
继 被 提出 。 目 前 ,直接 降低 管 电流 . 管 电 压 或 者 减少 扫 
描 角 向 采样 数目 是 实现 低 剂 量 CT 扫描 最 简单 且 最 有 效 
的 方法 。 但 是 ,降低 管 电流 会 导致 投影 数据 中 光子 噪声 
大 幅 增加 而 且 电 子 噪声 的 影响 也 会 更 加 突出 ;降低 管 电 
压 则 会 因 双 射线 的 穿 透 能 力 下 降 造 成 射 束 便 化 效应 ; 减 
少 扫描 角 向 采样 数目 使 得 投影 数据 不 满足 香农 采样 定 
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理 而 造成 重建 图 像 产 生 严重 的 条 形 伪 影 。 不 同 于 传统 
的 滤波 反 投影 算法 ,基于 低 剂量 CT 投影 数据 恢复 的 图 
像 重建 方法 对 CT 成 像 系统 进行 精确 建 模 ( 包 括 X 线 源 、 
探测 器 .电子 噪声 .成像 物体 等 ), 可 以 实现 优质 的 低 剂 
量 CT 图 像 重建 “人 ”。 其 目标 函数 是 根据 测量 数据 的 噪 
声 特性 来 构建 的 ,通常 包含 两 项 , 即 数 据 保 真 项 和 正则 
化 项 ,前 者 用 于 描述 投影 数据 的 统计 特性 ,后 者 用 于 修 
正解 的 特性 。 在 图 像 分 段 常 量 的 假设 下 ,全 变 分 (TV) 
正则 化 模型 * 其 在 去 除 噪 声 的 同时 可 以 保持 图 像 的 边 
缘 信息 。 然 而 ,许多 临床 CT 图 像 并 不 是 完全 满足 分 段 
常量 的 假设 ,在 剂量 特别 低 时 或 者 投影 角度 特别 少时 ， 
TV 重建 的 图 像 会 产生 阶梯 效应 和 块 状 伪 影 。 因 此 TV 
正则 化 在 基于 投影 数据 恢复 的 低 剂量 CT 重建 中 会 失去 
效力 。 为 了 解决 这 个 问题 ,本 文 提 出 一 个 基于 全 广义 变 分 
(TGV) 正 则 化 的 低 剂量 CT 重建 算法 ,其 可 以 在 投影 数 
据 不 满足 分 段 带 数 假设 的 前 提 下 去 除 噪 声 和 条 形 伪 影 。 
本 文 结构 安排 形式 如 下 :介绍 Anscombe 变换 公 
式 全 广义 变 分 概念 和 基于 全 广义 变 分 最 小 化 恢复 模型 
以 及 解 此 模型 的 优化 算法 ;对 Clock 体 模 和 Shepp- 
Logan 体 模 实验 结果 定性 和 定量 的 分 析 , 来 验证 所 提出 
方法 的 有 效 性 ;对 本 文 所 提出 的 方法 进行 总 结 和 讨论 。 


1 方法 
1.1 Anscombe € 3 

低 剂 量 CT 成 像 系统 中 ,探测 数据 可 以 描述 为 
I-I,exp(-u) ,其 中 工 为 X 射 线 入 射 光子 的 数目 ,1 为 出 
射 光子 的 数 日 并 目 服 从 Poisson 分 布 站 ,表示 经 过 对 数 
变换 后 的 投影 数据 , 即 Sinogram( 正 弦 图 )。 

服从 Poisson 分 布 的 数据 可 以 通过 Anscombe Z 
使 其 变 为 近似 服从 Gaussian 分 布 的 数据 。 假 设 x 服 从 
Poission 分 布 ,通过 如 下 的 Anscombe 变 换 : 


y(x)=2 xti (1) 


我 们 可 以 将 服从 Poission 分 布 的 x 变换 成 近似 服从 
JT 2:8 MI] Gaussian 4) ff^", 

假设 投影 数据 每 个 像素 间 是 互相 独立 的 ,并 且 经 过 
Anscombe 变 换 后 都 服从 方差 为 1 的 Gaussion 分 布 , 令 
fe se) 为 待 估计 投影 数据 的 理想 值 , 
u — (uu, uu). 表示 投影 数据 经 过 Anscombe 变 化 后 
近似 满足 Gaussian 分 布 的 数据 ,其 中 表示 转 置 运 算 。 
1.2 全 广义 变 分 (TGV) 

TGV 正则 化 模型 由 Berdies 557" 7 提出 并 用 于 网 像 
去 品 。 给 定 图 像 fk 阶 TGV 的 定义 为 : 

TGVf)- sup [,u div' vdx|v e C(Q, Sym (R^) , 

Ilvll 5a. ]ldiv vll, 5) v) 


其 中 ,1=0,…,k-lLkeN & TGV 的 阶 数 ， 


æ 


a= (a, 0, ts Ori) ÆTGV AJI E, Sym (R^) 是 
阶 对 称 张 量 空 间 , ! 阶 散 度 算 子 定义 是 : 
| 1! OV,, 
divv-2,,,—— NEM, (3) 
U* Ox 
其 中 , M,2(eN'IE2 mk) ,大 阶 张 量 对 称 张 量 
的 无 穷 范 数 定义 为 : 


Ivl], = sup | Y Ever (4) 


M k=1 Hf, TGV 和 TV 是 等 价 , 即 可 以 表示 为 : 
TGV*=a6TV 。 在 本 文中 ,我 们 只 考虑 二 阶 TGV ,定义 
为 : 

TGV;(f)- supt[,u div'vdx|v e CQ, S^) ， 

| vll, <alldiv vll, Sa; 

其 中 , S$ 是 dx4q 维 矩 阵 空间 , a= (a6。，Q) 是 
正 数 ,进一步 ,对 称 和 矩阵 的 一 阶 和 二 阶 散 度 可 以 由 下 式 
表示 : 

(div v), Ya dicis > E D A (6) 


j i-l 


类 似 的 ,(5) 式 中 的 无 穷 范 数 定义 为 : 
-sp (Sof — m 


| (8) 


其 中 ,文献 ”给 出 了 与 式 (5) 互 相等 价 的 二 阶 TGV 
的 定义 : 

TGVAf)- mina, fal Vf- o|dx +a, fleo) ax (9) 

其 中 , e(o) =A Vo +Vo") H Radon dE, a, fl a, 
是 TGV 的 权重 ,来 平衡 (5) 式 的 前 一 项 的 一 阶 导数 和 后 一 
项 的 二 阶 导数 。 由 于 TGV 用 高 阶 导数 信息 描述 图 像 的 特 
性 ,TGV 不 需要 满足 投影 数据 是 分 段 常 数 的 假设 >。 
1.3 基于 投影 数据 TGV 最 小 化 的 CT 重建 方法 

投影 数据 恢复 的 一 般 模 型 他 可 以 表示 为 : 

f- arg min Ju -fÈ +R) (10) 


其 中 ,第 1 项 为 数据 保 真 项 ,第 2 项 为 正则 化 项 , 
4>0 是 正则 化 参数 。 

本 文 用 TGV 正则 化 代替 式 (10) 中 的 正则 化 项 
R( 了 ) 得 到 基于 TGV 最 小 化 的 理想 数据 恢复 模型 为 ; 

f- arg min dlu -fÈ -2TGY (7) (11) 


其 中 ,是 正则 化 参数 。 

公式 (11) 的 离散 形式 为 : 

min, Hlu-f IEAA- ol Ale), — Q2 

其 中 , FÆR" W=R®" 以 及 微分 算 子 div , e, V 
可 以 由 限 差分 算 子 近 似 得 到 ”1。 


(5) 


a se OV, 
|div v|| =sup 3$ ! (x) 
i| 71 Xj 


J 
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根据 对 偶 原 理 2 ,我 们 可 以 将 (12) 转 化 为 如 下 的 鞍 
点 问题 : 
min max EIL -FIÈ +D Vf-œ@)+q"el@) (13) 


feFweW peP.qeQ 


其 中 , p 和 g 为 对 偶 变 量 , 并 且 忆 和 2 定义 为 : 


P-ipeR"""|p|| x53 (14) 
Q-iq e R""[lal, S5) (15) 
proj (P) fl proj,(q) 的 定义 如 下 : 
prof (B) - s [BIB) (16) 
proj(ã)= RUI (17) 
Proxl( 刀 是 邻近 算 子 ,其 由 下 式 表示 : 

uf JE 
prox) aremin 2h + | X Z (18) 


由 上 述 方法 求 出 (11) 式 的 解 并 对 其 进行 Anscombe 
逆 变 换 , 最 后 由 滤波 反 投 影 算法 进行 重建 。 
下 面 我 们 给 出 基于 TGV 正则 化 的 低 剂量 CT 重建 
方法 的 计算 步骤 : 
步 又 1 对 低 剂 量 CT 投 影 数据 进行 Anscombe 变 换 ， 
将 服从 Poission 分布 的 投影 数据 转化 为 近似 Gaussian 
分 布 的 数据 。 
步骤 2 基于 全 广义 变 分 最 小 化 模型 ,我 们 使 用 
Chambolle-Pock 算 法 2 求解 ,从 而 达到 对 Anscombe 变 
换 后 投影 数据 去 噪 的 目的 。 
() 初 始 化 参数 : Bos BiT p3 
(2) 初 始 化 变量 : p, d,o o f^. f 
(3) 对 数据 进行 迭代 恢复 ,迭代 形式 如 下 : 
For n=0,1,°…,N 
p" -proj p" *p(Vf* - a) ; 
q''-projj(q' * pe(a/)) ， 
y 
f" 2 proxi(f" * cdiv,p"'), 
JESAS m 


old _ n 
eo =w 


n+ 


œ" ' =œ" +(p"'+div,q"') ， 
à -2g" - a", 
End 
步骤 3 对 步 又 2 得 到 的 数据 其 进行 Anscombe 逆 变 
换 , 得 到 去 噪 后 的 投影 数据 。 
步骤 4 对 步骤 3 得 到 的 投影 数据 进行 EBP 重建 得 
到 CT 图 像 。 


2 实验 与 分 析 
为 了 验证 全 广义 变 分 最 小 化 方法 在 低 剂量 CT 图 


像 重 建 中 的 有 效 性 ,本 文通 过 Clock 和 Shepp-Logan 数 
值 体 模 进 行 定量 和 定性 分 析 。 并 且 与 斜 波 滤波 反 投影 
法 . 汉 宁 滤波 反 投 影 算 法 和 非 单调 全 变 分 最 小 化 算法 加 
进行 比较 。CT 成 像 几 何 采 用 弧 形 探测 器 的 扇形 束 ,其 
中 射线 源 到 旋转 中 心 和 探测 器 的 距离 分 别 为 570 mm 
和 1040 mm, 旋 转角 在 [0, 2r] 间 的 采样 值 为 1160, 探 测 
器 个 数 为 672。 体 模 大 小 设 定 为 312x512, 对 于 Clock 
体 模 ,入 射 光 子 数 为 $.0x104; 相 应 的 Shepp-Logan 体 模 
的 入射 光子 数 为 1.0x10;。 
2.1 Clock 体 模 实验 

Clock 体 模 实验 结果 如 图 1 所 示 。 其 中 ,图 1A 为 真 
实 的 体 模 图 像 ; 图 1B 为 斜 波 滤波 反 投 影 算法 重建 的 结 
果 ; 图 1C 为 汉 宁 滤波 反 投影 算法 重建 的 结果 ;图 1D 为 
韭 单调 全 变 分 最 小 化 (NTVM) 算 法 重建 的 结果 ;图 1E 
为 本 文 提出 的 全 广义 变 分 最 小 化 方法 重建 的 结果 。 为 
了 进一步 展现 本 方法 重建 效果 ,图 2 给 出 了 对 应 图 1 的 
局 部 放大 图 。 由 图 2 可 以 清晰 的 看 出 本 文 方法 对 噪声 
和 伪 影 的 抑制 效果 。 

通过 计算 重建 图 像 的 信 噪 比 和 均 方 误差 可 以 定量 
地 分 析 不 同方 法 的 差异 。 信 了 噪 比 计算 公式 为 : 


SNR=10 oed 2 (F..()-F..) | (19) 
AF.) - FA) 
均 方 误差 计算 公式 为 : 
YF. )-F..0) | vd 
2 (F1) 

其 中 , FL) 表示 像素 点 i 通过 算法 得 到 的 重建 图 
像 的 灰 度 值 , ,表示 重建 图 像 已 . 中 所 有 像素 点 灰 度 
的 平均 值 ; Fueli) 表示 真实 图 像 像 素 点 i 处 的 灰 度 值 。 

HX 1 的 性 噪 比 和 均 方 误差 分 析 可 以 看 出 ,相对 于 
滤波 反 投 影 和 非 单调 最 小 化 方法 ,本 文 提 出 的 重建 图 像 
的 方法 具有 最 高 的 信 噪 比 和 最 小 的 均 方 误差 。 与 
NTVM 算 法 相 比 ,TGV 方 法 信 噪 比 提高 了 2.619% ,而 相 
应 的 均 方 误差 减少 了 13.04%。 

为 了 进一步 验证 本 文 提出 方法 的 优越 性 ,我们 分 别 
给 出 了 4 种 方法 重建 结果 的 线 剖面 图 (图 3) ,其 中 横 坐 
标 为 图 像 像素 位 置 , 横 坐 标 为 图 像 的 衰减 值 。 由 图 3 可 
以 看 出 本 文 提 出 的 重建 图 像 的 方法 与 其 他 3 种 方法 相 
比 ,能 够 较 好 的 接近 真实 图 像 的 线 剖面 。 

2.2 Shepp-Logan 体 模 实 验 

Shepp-Logan 体 模 实验 结果 如 图 4 所 示 。 其 中 ,图 
4A 为 真实 的 体 模 图 像 ;图 4B 为 通过 和 斜 波 滤波 反 投影 
法 重建 的 结果 ;图 4C 为 汉 宁 滤波 反 投影 算法 重建 的 结 
果 ; 图 4D 为 非 单调 全 变 分 最 小 化 算法 重建 的 结果 ;图 
4E 为 本 文 提 出 的 全 广义 变 分 方法 重建 的 结果 。 网 5 给 
出 了 低 剂量 CT 重建 结果 的 一 个 局 部 的 放大 图 。 


NMSE = | 


C 
图 1 Clock 体 模 图 像 和 不 同方 法 重建 的 图 像 


© 


Fig.1 Clock phantom and image reconstructed by different methods. A: Clock phantom; B: Image reconstructed 
by FBP algorithm with ramp filter; C: Image reconstructed by FBP algorithm with hanning filter; D: Image 
reconstructed by NTVM method; E: Image reconstructed by TGV method. 


图 2 Clock 体 膜 图 像 局 部 放大 图 

Fig.2 Zoomed- in views of the Clock phantom and images 
reconstructed by different methods. Clock phantom (first 
column); Image reconstructed by FBP algorithm with ramp filter 
(second column); Image reconstructed by FBP algorithm with 
hanning filter (third column); Image reconstructed by NTVM 
method (fourth column); Image reconstructed by TGV method 
(fifth column). 


X1 Clock 体 膜 图 像 的 信 噪 比 和 均 方 误差 
Tab.1 SNR and NMSE of the Clock phantom 


Methods FBP (ramp) FBP (hann) NTVM TGV 
SNR (dB) 17.7521 22.0859 23.422] . 24.0352 


NMSE 0.0086 0.0031 0.0023 0.0020 


由 表 2 的 信 噪 比 和 均 方 误差 分 析 可 以 看 出 ,与 其 他 
3 种 方法 相 比 ,本 文 提出 的 全 广义 变 分 方法 在 提高 图 像 
重建 的 信 噪 比 ,降低 重建 图 像 均 方 误差 方面 均 有 上 佳 表 
现 。 与 NTVM 算 法 相 比 ,TGV 方法 信 噪 比 提高 了 
1.03% ,而 相应 的 均 方 误差 减少 了 8%。 

此 外 ,为 了 更 直观 的 验证 重建 图 像 的 质量 ,我 们 画 
出 了 重建 图 像 的 线 前 面 图 (图 6) ,其 中 横 坐 标 为 图 像 像 
素 位 置 , 横 坐 标 为 图 像 的 衰减 值 。 由 线 剖 面 图 可 以 看 出 
本 文 提出 的 方法 重建 的 图 像 与 真实 图 像 的 线 剖 面 有 更 
高 的 吻合 度 。 


3 讨论 

基于 TV 最 小 化 的 CT 图 像 重 建 方法 在 抑制 噪声 和 
伪 影 方面 取得 了 很 好 的 效果 。 由 于 TV 是 基于 图 像 分 段 
常量 的 假设 , 当 品 声 水 平 非常 高 时 ,会 产生 严重 的 阶梯 
效应 。 针 对 TV 方法 重建 的 图 像 会 产生 阶梯 效应 的 问 
题 ,本 文 提出 了 一 个 基于 TGV 正 则 化 的 低 剂 量 CT 重 建 
算法 ,该 方法 不 需要 在 投影 数据 满足 分 段 常数 假设 。 由 
于 TGV 利 用 了 高 阶 导 数 信息 ,在 不 满足 分 段 常量 的 假 
设 下 ,本 文 方法 可 以 有 效 去 除 噪声 的 同时 ,可 以 有 效 地 
抑制 噪声 与 条 形 伪 影 。 数 值 仿真 实验 结果 表明 本 文 提 
出 的 方法 可 以 有 效 地 抑制 低 剂 量 CT 图 像 中 的 噪声 和 伪 
影 ,消除 TV 正常 产生 的 阶梯 效应 ,从 而 改善 低 剂 量 CT 
图 像 重 建 的 质量 。 
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图 3 Clock 体 膜 图 像 的 线 剖 面 图 
Fig.3 Profiles of the clock phantom. 


图 4 Shepp-Logan 体 模 图 像 和 不 同方 法 重建 的 图 像 


Fig.4 Shepp-Logan phantom and images reconstructed by different methods. A: Clock phantom; B: 


Image reconstructed by FBP algorithm with ramp filter; C: Image reconstructed by FBP algorithm 
with hanning filter; D: Image reconstructed by NTVM method; E: Image reconstructed by TGV 
method. 


在 Clock 体 膜 实验 中 ,由 图 1 可 以 看 出 和 斜 波 滤波 反  ” 波 反 投影 和 非 单调 全 变 分 最 小 化 方法 ,我 们 提出 的 全 广 
投影 和 汉 宁 滤波 反 投影 算法 重建 的 图 像 含 有 大 量 的 噪 义 变 分 最 小 化 方法 在 抑制 噪声 与 条 形 伪 影 和 保持 边缘 
声 和 条 形 伪 影 ,图 像 质 量 严重 退化 。 非 单调 全 变 分 算法 ”方面 都 有 良好 表现 ,可 以 更 好 的 保持 图 像 一 致 性 。 进 一 
重建 的 图 像 虽 然 在 一 定 程度 上 减少 了 CT 图 像 的 条 形 伪 ”  ” 步 ,由 体 膜 图 像 重 建 结果 的 放大 图 可 以 明显 的 看 出 ,本 
影 和 噪声 ,但 由 于 TV 正则 化 的 缺陷 ,部 分 伪 影 和 噪声 仍 ”” 文 提出 的 全 广义 变 分 方法 重建 的 图 像 对 噪声 和 伪 影 
然 存在 。 从 视觉 效果 评价 来 看 ,我们 提出 的 方法 在 去 除 。 ”更 好 的 抑制 作用 。 此 外 ,Clock 和 Shepp-Logan 体 膜 实 
图 像 的 噪声 和 条 形 伪 影 方面 具有 更 佳 的 表现 。 在 。”” 验 的 定量 实验 结果 进一步 表明 ,本 文 提出 的 方法 比 非 单 
Shepp-Logan 体 膜 实验 中 ,由 图 4 我 们 可 以 看 到 , 相 比 滤 调 全 变 分 算法 更 具 优 越 性 。 
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图 5 Shepp-Logan 体 膜 图 像 局 部 放大 图 
Fig.5 Zoomed-in views of the Shepp-Logan phantom and images reconstructed by different methods. A: 
Clock phantom; B: Image reconstructed by FBP algorithm with ramp filter; C: Image reconstructed by FBP 
algorithm with hanning filter; D: Image reconstructed by NTVM method; E: Image reconstructed by TGV 
method (fifth column). 


5&2 Shepp-Logan 体 膜 图 像 的 信 噪 比 和 均 方 误差 
Tab.2 SNR and NMSE of the Shepp-Logan phantom 


由 于 本 文 提出 的 方法 需要 根据 具体 的 重建 目标 人 
工 选择 正则 化 参数 ,大 大 限制 了 其 更 为 广泛 的 应 用 。 从 


Methods FBP (amp) FBP(hann) NIVM TGV 国内 外 的 研究 现状 来 看 ,目前 还 没 一 种 有 效 的 参数 选择 
SNR (dB) 19.3790 22.2822 23.1767 23.4181 方式 ,使 得 正则 化 的 自 适 应 参数 调整 达到 最 住 状态 。 
NMSE 0.0058 0.0030 0.0025 0.0023 此 ,如 何 自 适应 选择 正则 化 化 参数 是 我 们 今后 研究 的 一 
个 重点 课题 。 
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图 6 Shepp-Logan 体 膜 图 像 的 线 剖 面 图 
Fig.6 Profiles of the Shepp-Logan phantom. 
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